Large Broadening of the Superconducting Transition by Fluctuations in a 3D Metal at 

High Magnetic Fields: The MgB 2 case. 
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It is shown that the transition to the low temperature superconducting state in a 3D metal 
at high magnetic field is smeared dramatically by thermal fluctuation of the superconducting order 
parameter. The resulting superconducting-to-normal crossover occurs in a vortex liquid state which 
is extended well below the mean- field H C 2. Application to MgB2 yields good quantitative agreement 
with recently reported data of dHvA oscillation in the superconducting state. 



It is well known that the transition from the normal 
to the superconducting (SC) state in typc-II 3D super- 
conductors in the absence of external magnetic field is 
a sharp, second-order phase transition, with a vanishing 
order parameter at the transition temperature T c contin- 
uously growing with the decreasing temperature below 
T c . Fluctuations effect can smear the transition sig- 
nificantly in high T r and low-dimensional superconduc- 
tors 0|, where the phase space accessible for the fluc- 
tuations is dramatically enhanced. The influence of an 
external magnetic field is similar to an effective reduc- 
tion of dimensionality 0, resulting in a significant smear- 
ing of the transition even at very low temperatures. Such 
strong smearing effects have been observed in varies quasi 
2D low T c superconductors at high magnetic fields 

In the present paper we show theoretically, and con- 
firm by comparison with very recent de-Haas van-Alphcn 
(dHvA) oscillation data in the SC state @, that the 
smearing of the SC transition by fluctuations in a con- 
ventional 3D type-II superconductor, such as MgB2, at 
high magnetic fields, is surprisingly strong, comparable 
in magnitude to that in 2D superconductors. This con- 
clusion is reached by generalizing the Bragg-chain model 
of the 2D vortex liquid state at high perpendicular mag- 
netic field 01 to an array of strongly coupled 2D SC layers. 
Within this model it is found that, similar to the situa- 
tion in a single 2D SC layer, the vortex lattice melting 
point in a 3D superconductor at low temperature T is lo- 
cated well below the mean-field (MF) upper critical field 
H C 2 (T), so that in a broad field range above the melting 
point the corresponding system of fluctuations is equiva- 
lent to a ID array of SC quantum dots at zero magnetic 
field @. 

Our starting point is the microscopic BCS Hamiltonian 
for electrons in a layered 3D metal, interacting via an ef- 
fective two-body attractive potential, under the influence 
of a strong static magnetic field. We assume, for sim- 
plicity, that the magnetic field direction (along z-axis) is 
perpendicular to the layers situated in (x, y)-plane. Writ- 
ing down the functional integral expression for the par- 
tition function of this system, the electronic field can be 
eliminated by introducing bosonic complex field A (r) ( 



by means of the Hubbard-Stratonovich transformations), 
which describes all possible realizations of Cooper-pairs 
condensates 0,0- Expansion of the resulting free en- 
ergy functional, Fq [A (r)], in the order parameter up to 
the quartic term is a good approximation for magnetic 
fields around mean field iT C 2- 

In the lowest Landau level approximation, which is 
valid at high magnetic fields and low temperatures [lflj . 
the most general form of the order parameter A (r) 
is a coherent superposition of Landau wave functions, 
(j> q (x,y) = exp iqx - (y/a H + qa H /2Y 



A (x, y,z) = ^2 c q ( z ) ^q ( x > v) 



(1) 



where c q (z) are arbitrary complex functions of the co- 
ordinate z , defined on the quasi continuous lattice 

1 = iB^i ( n + m/y/N) ,n,m= -y/N/2 + 1, . . . , VN/2, 
which determines the projections of the orbital guiding 
centers on the y axis. Here N is the total number of flux 
lines threading the SC sample, and ajj = y ch/eH is the 
magnetic length. In the case when all the coefficients c q 
arc different from zero there is one to one correspondence 
between all N guiding centers and their projections on 
the y axis. 

As discussed in detail in Ref.|lll|. such a quasi contin- 
uous configuration costs a large fraction of the total SC 
condensation energy, and may be therefore omitted at 
magnetic fields well below MF H C 2- In this region it is 
sufficient to use in Eq. (JTJ a discrete subset of \/~N Landau 
orbitals <f> qn (x,y), with q„ 



: n, n = 0,±I,±2, . 



u q n v~>*/_> a x a H ' 

Including the omitted high energy modes would always 
increase the fluctuation effect. Thus, near and above the 
MF transition the fluctuation effect, calculated in the 
discrete chain approximation, should always underesti- 
mate the observed effect. On the other hand, an upper 
bound on the fluctuation effect can be obtained from the 
2D model ^3], where the discrete chain representation 
yields reasonably good agreement with exact numerical 
simulation [l^ . 



Writing c qn (z) 



(*)l 



the phase 



if n (z) determines the relative lateral position x n = 
—ifn/qn of the n-th Landau orbital within a single 2D 



2 



layer. It can be readily shown that, by selecting the x 
axis along the principal crystallographic axis of the 2D 
Abrikosov vortex lattice, the variables £„ = tp n — f n -i 
describe the lateral positions of the most easily sliding 
vortex chains in the vortex lattice, which are generated 
mainly by interference between two neighboring Landau 
orbitals [nj. 

The corresponding free energy functional, projected on 
the subspace of the lowest Landau level, can be written in 
a local anisotropic Ginzburg-Landau (GL) form 0,E3 : 



Fi 



GL 



d 3 r 



-a\A{r) 



f |A(r)| 4 + 7 



3 A (r) 



dz 



(2) 

where v = ira^d is the volume of a single vortex per 
layer and the effective GL coefficients a, ft, and 7, can 
be expressed in terms of the microscopic normal elec- 
tron parameters (see below). It is convenient to divide 
the length L z of the sample perpendicular to the lay- 
ers into N z segments of length d , L z = N z d , where 
d is the intcrlayer distance, and define a discrete set 
c n , C = I e^-e = c„(z c ) , C = 1,2, ...iV, , with the 
periodic boundary conditions: c u _n z +i = c n ,i- 
Using Eq. (|TJ the partition function, 

Z = Z^ = J DA (r) DA* (r) exp {-F GL [A (r)] /k B T} 

where Z c h is evaluated as a multiple integral: 

n„, c / KcMKd / di PnX e- F <^/(VNk B T) ) with 



(a x / 

h 2 3/2 



^ 7 = = -«EKd 2 + */EKc +1 -^ c l 



E a 



2_ L(ll 2 . 



xe i(^+s,c+^+p,c- ( P™,c- i ^+s+P,<;) 



(3) 



and r\ = d 2 ~f. In this Landau orbital representation of 
the GL functional, Eq. , the off-diagonal elements con- 
stitute a rapidly convergent series in the small expansion 
parameter A = e~^' ax ' w 0.066, where all terms of or- 
der higher than the second may be neglected. 

For the 3D electronic band structure under study here 
the characteristic interlayer Josephson tunneling ampli- 
tude 77 is much larger than the minimal intra-layer shear 
stiffness, 4A 2 ~ 10~ 2 , characterizing the principal axis, 
x, in the vortex lattice. Under this condition the low en- 
ergy fluctuations of A (r) correspond to collectively slid- 
ing chains of vortices in different layers, such that the cor- 
responding fluctuating magnetic flux lines remain nearly 
parallel to each other (and so to the external magnetic 
field). Similar to the situation in a 2D superconductor 
, these low- lying shear fluctuations determine the vor- 
tex lattice melting point to be well below H C 2. 

For magnetic fields H above this melting point the 
second order terms in A, which oscillate as functions of 



the phases <p n ,c, are averaged to zero by the integrations 
over <£>n,c, an d the effective free energy functional Fql 
can be written in the very simple, independent vortex 
chain form F GL = VNJ2 n foL > witn 



+ V \Cn,X+l ~ C «,C| 2 } i 



-^-Ic r\ 4 

2 3/2 l ClI <Cl 



(4) 



which may be considered as an effective GL energy func- 
tional for a single vortex line. In this expression we have 
also neglected the first order terms in A, since they effec- 
tively yield a small additive correction to (3, so that their 
influence on the critical behavior is unimportant |l3j . 

The calculation of the corresponding partition function 
is rather straightforward. The phase variables, <p n £, can 
be readily integrated out. The resulting expression can 
be written as a functional integral over the squared am- 
plitude y^ = (Pa x /2kBTy/n) 1 ^ 2 \c^\ 2 , with an effective 
GL functional, incorporating phase fluctuations: 

fGL.eff = ^b^E ( - V2xy c + iy 2 

C 



In J (2«VZ/cyCH 



2ny c 
(5) 



n 



Here and below, for the sake of notation simplicity, we 
drop the vortex chain indices. 

We can then estimate the partition function in two lim- 
iting situations, for weak (re <C 1) and strong (re 3> 1) 
interlayer coupling. In the former limit the integrals 
over amplitudes can be calculated explicitly with the well 
known result Q,E3- It is convenient to define the par- 
tition function per single vortex per layer, Z v = Z 1 ^ 
(TV = NN Z ), which is given by: 



In erfc 



(6) 



where erfc(x) = -|= e v dy. For strong interlayer 
coupling the result is obtained by using the steepest de- 
scend integration, leading to 

Z 11 

In = V2xy - -vl - - In (4ttk) 



,V4 



•In- 



1) 



1/2 



2V2, 



2 3 / 4 



(7) 



where y Q = (x + Vx 2 + l) /y/2. In Eas. pTj) Z and Z 
are constants (i.e. independent of both x and re), and so 
are thcrmodynamically unimportant. 

Thus, in the liquid state above the vortex lattice melt- 
ing point one can readily derive simple limiting expres- 
sions for the spatially averaged mean square order param- 
eter, (|A| 2 \ = / d 3 r(\A(T*)\ 2 ) /(Afv), in the general 
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FIG. 1: The spatially averaged mean-squared order param- 
eter, measured in units of aoi , as a function of the scaling 
parameter x = a/ ^2/3p a ksT (see text ) around the mean- 
field transition point x — in the weak and strong coupling 
limits (Solid lines). Dashed lines represent the corresponding 
results calculated from a simple interpolation formula. The 
doted straight line represents the result of mean-field theory. 
Inset: The derivative with respect to x is plotted to emphasize 
the crossover region. 
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where evo/ 



2/3/3 a - 



The function $0(2;; k) 



a i$ (x; k) 
(8) 

: for 

ox 



K = and k = 100 is shown in Fig. (__). With increas- 
ing interlayer coupling the mean-square order parame- 
ter at the MF transition point decreases rapidly within 
the small interval < k < 1, then saturating at a non- 
vanishing value in the entire strong coupling region k > 1. 

at the MF tran- 



The non- vanishing value of (JA| 
sition point reflects smearing of the phase transition to 
the SC state due to fluctuations. This effect characterizes 
the GL theory of a 3D superconductor at high magnetic 
fields, since above the vortex lattice melting point it can 
be mapped to the GL theory of a ID superconductor 
at zero magnetic field, where a genuine phase transition 
is absent 1 1 5| . In the zero coupling limit, rj — > , the 
effective GL free energy functional in a given SC layer, 
appearing within the curly brackets in Eq. Q , is equiva- 
lent to the GL energy functional for a 0D superconductor 
near the SC transition at vanishing magnetic field 0: 



f (0) 

J GL 



(Vi 2 ) 



k B T c 



i-O 



2 . 0-106 . ,4 



Here T c0 = T C (H -»■ 0), t = T/T c0 , and S = 
[N (0)uA: B r c o] _1 ( with N (0)- the electronic density of 



states per unit volume at the Fermi energy, and v- the 
volume of the small SC particle), is the quantum size 
parameter of the SC grain. The corresponding MF free 
energy is f§l (T,H->0) = - ( $£) k B T c0 (1 - tf . 

For the isolated SC layer at high magnetic field and low 
temperature T<T C (H — > 0) it was found that _jfl , a = 

i^M^^itM 1 -^™* 11 ^^/^ (T-0) , 

and P f» 7ft~lrg ■ Here ui c = eH/m*c is the cyclotron 
frequency, with m* the in-plane effective cyclotron mass, 
and Ep the Fermi energy. The corresponding MF 



-m ff L {T,H) = /g_(0,0)(l 

2 n tt / tt trn\i2 



fgl(Q,0)(l-T/T c0 r[l-H/H c2 (T)] 



free energy is (T —> 0,H) w — (^J £ F (1 - ft) 2 . 
For general T , and H values near the MF transi- 
tion line, H c2 (T) = = H c2 (0) (1 - T/T c0 ), 
the corresponding MF free energy has the well 

t-hf = 
so that by 

equating the coefficient Jq L (0, 0) in the two limit- 
ing regions of the phase boundary, i.e. /_?__ (0, 0) = 
(f|)Mo = (t^)^, one finds: 8 = 
with f3 a « 1. 

Using the well known BCS expression for the zero tem- 
perature coherence length, £ (0) = 0.18hvp /k B T c0 , we 
find that 



28 



1 



(9) 



W(0), 

This enables us to estimate the effective spatial size of 
the SC grain in the zero field limit. Since our original 
model system consists of a 2D SC layer, the volume v = 



______ 

4tt Ef : 



so that S 



electron gas, N (0) = m*/2-Kh 2 

35 (^a^) (fcjbj) ' Comparing this expression to Eq.© , 
we find that the radius of the effective SC grain a ~ £ (0) 



which is approximately equal to Qh 



2(0)" 



the smallest 



length scale in a 2D SC condensate in a magnetic field 
H w ff c2 (0). For MgB 2 we have: T c0 w 40if , E F w 
12500K, so that <5 = 0. 23 . 

These results show that a single SC layer in a high 
magnetic field above the vortex lattice melting point is 
equivalent to a 0D superconductor at zero magnetic field 
|8j. The coupled layer model is therefore equivalent to a 
ID Josephson array of small SC grains, discussed, e.g. in 
or in _________ . For the relatively large effective quantum 

size parameter 5 = 0. 23 , estimated above, the numerical 
simulations reported in Ref.H show that the reduction 
of the transition width by the interlayer coupling is not 
very significant. 

To make this feature useful in our analysis of experi- 
mental data we have approximated our result by means 
of an interpolation function, suggested by Ito et al.jjjj to 
fit their experimental dHvA data for a quasi 2D organic 
superconductor, namely: 



interp 



(x; k) = x + 



(10) 



4 



where the fitting parameter v K depends only on the inter- 
layer coupling. The equivalence of Eq. l|10|l to the inter- 
polation formula presented in Ref.[j| is made clear if we 

note from Eq.© that x = = ^ (l - h^t)) , 

where Ao is the SC gap parameter at T = 0, and H = 0. 

The best fitting function ^o,i n terp{x; k) is represented 
in Fig. Q by the dashed lines for the limiting cases k = 
and k — > 00 . At zero coupling the parameter v K can 
be obtained by comparing &o,interp(x', 0) with &o(x; 0) = 
2 (x + exp(— x 2 ) / y/ir erk (— x)) , to yield u K= o = 2/y/n « 
1.13. In the strong coupling limit the best fit is obtained 
for v K= ioo = -51. Thus, the entire range of the actual 
smearing parameter ai = v K aoi is obtained with v K = 
.51 - 1.13. 

The parameter aj = v K ^J -^pjf- > which controls the 
smearing of the phase transition by the fluctuations, is 
related to the parameter ap introduced in Ref. [|| (de- 
noted there by a ) by: qj? = 2a// A§ , where Ao can be 
identified with Ae of Ref. @|. Using the above value 
of f3, obtained from the 2D electron gas model, we find 
that: 

aj m 0.35v K huj c y // Epkp,T (11) 

which reflects the smearing effect due to the magnetic 
field, the suppression of this smearing by the interlayer 
coupling parameter ft, and the thermal smearing effect. 

For the parameters characterizing the MgB2 data, re- 
ported by Fletcher et al. |(|, i.e. with m* m 0.3m e 
, T = 0.32A , A = 200A, F = 2930T , so that 
Ep = hui c F/H = 12464A, our theoretical estimate is 
aj m A79K 2 , or ap f» 0.024. This result is similar to, 
but somewhat smaller than the best fitting value of ap 
obtained in Ref.0. 

This is quite a reasonable result since, as noted above, 
one should regard the calculated width as a lower bound 
on the smearing of the SC phase transition by all types 
of thermal fluctuations of the SC order parameter. Fur- 
thermore, as clearly seen in Fig.JTJ, the deviation of the 
fluctuation effect predicted in our 3D model from the best 
fit to the experimentally observed data is rather small in 
the magnetic field region bellow MF H C 2, and becomes 
more significant at fields well above the MF transition. 
This behavior is reasonably explained by the omission of 
many fluctuation degrees of freedom in our vortex chain 
model ^3 1 which is expected to weaken progressively the 
fluctuation effect with respect to the actually observed 
one as the magnetic field increases above MF 

It should be also emphasized that the smearing param- 
eter cti, given in Ea. (|llf) . is independent of the SC gap 



parameter Ao, an important parameter determining the 
damping of the dHvA oscillations in the SC state |j| ■ In 
the fitting procedure employed in Rcf.Q Ao was found 
to disagree with the SC gap parameter A„ derived by 
other methods. Our estimate of the SC fluctuations ef- 
fect from the damping of the dHvA oscillation in the SC 
state is thus unaffected by such a disagreement. 

In conclusion, it was shown here that the transition to 
the low temperature SC state in a pure, extremely type- 
II, 3D superconductor is smeared dramatically by the 
magnetic field, comparable in magnitude to the broad- 
ening of the transition observed in quasi 2D supercon- 
ductors. The dimensionality reduction by the magnetic 
field, responsible for this smearing, is shown to take place 
in a broad field range above the vortex lattice melting 
point, where the system of SC fluctuations is equivalent 
to a ID array of SC quantum dots. The theoretically 
predicted width of the transition region is found to be 
in good agreement with experimental data of the dHvA 
effect in the SC state of MgB 2 . 
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